Clinical and epidemiological characteristics of Madariaga and Venezuelan equine encephalitis virus infections.

Madariaga virus (MADV) and Venezuelan equine encephalitis virus (VEEV) are emerging arboviruses affecting rural and remote areas of Latin America. However, there are limited clinical and epidemiological reports available, and outbreaks are occurring at an increasing frequency. We addressed this gap by analyzing all the available clinical and epidemiological data of MADV and VEEV infections recorded since 1961 in Panama. A total of 168 of human alphavirus encephalitis cases were detected in Panama from 1961 to 2023. Here we describe the clinical signs and symptoms and epidemiological characteristics of these cases, and also explored signs and symptoms as potential predictors of encephalitic alphavirus infection when compared to those of other arbovirus infections occurring in the region. Our results highlight the challenges clinical diagnosis of alphavirus disease in endemic regions with overlapping circulation of multiple arboviruses.

Tracking arboviruses, their transmission vectors and potential hosts by nanopore sequencing of mosquitoes.

The risk to human health from mosquito-borne viruses such as dengue, chikungunya and yellow fever is increasing due to increased human expansion, deforestation and climate change. To anticipate and predict the spread and transmission of mosquito-borne viruses, a better understanding of the transmission cycle in mosquito populations is needed. We present a pathogen-agnostic combined sequencing protocol for identifying vectors, viral pathogens and their hosts or reservoirs using portable Oxford Nanopore sequencing. Using mosquitoes collected in São Paulo, Brazil, we extracted RNA for virus identification and DNA for blood meal and mosquito identification. Mosquitoes and blood meals were identified by comparing cytochrome c oxidase I (COI) sequences against a curated Barcode of Life Data System (BOLD). Viruses were identified using the SMART-9N protocol, which allows amplified DNA to be prepared with native barcoding for nanopore sequencing. Kraken 2 was employed to detect viral pathogens and Minimap2 and BOLD identified the contents of the blood meal. Due to the high similarity of some species, mosquito identification was conducted using blast after generation of consensus COI sequences using RACON polishing. This protocol can simultaneously uncover viral diversity, mosquito species and mosquito feeding habits. It also has the potential to increase understanding of mosquito genetic diversity and transmission dynamics of zoonotic mosquito-borne viruses.

Chikungunya: a decade of burden in the Americas.

In the Americas, one decade following its emergence in 2013, chikungunya virus (CHIKV) continues to spread and cause epidemics across the region. To date, 3.7 million suspected and laboratory-confirmed chikungunya cases have been reported in 50 countries or territories in the Americas. Here, we outline the current status and epidemiological aspects of chikungunya in the Americas and discuss prospects for future research and public health strategies to combat CHIKV in the region.

Real-time RT-PCR for Venezuelan equine encephalitis complex, Madariaga, and Eastern equine encephalitis viruses: application in human and mosquito public health surveillance in Panama.

Eastern equine encephalitis virus (EEEV), Madariaga virus (MADV), and Venezuelan equine encephalitis virus complex (VEEV) are New World alphaviruses transmitted by mosquitoes. They cause febrile and sometimes severe neurological diseases in human and equine hosts. Detecting them during the acute phase is hindered by non-specific symptoms and limited diagnostic tools. We designed and clinically assessed real-time reverse transcription polymerase chain reaction assays (rRT-PCRs) for VEEV complex, MADV, and EEEV using whole-genome sequences. Validation involved 15 retrospective serum samples from 2015 to 2017 outbreaks, 150 mosquito pools from 2015, and 118 prospective samples from 2021 to 2022 surveillance in Panama. The rRT-PCRs detected VEEV complex RNA in 10 samples (66.7%) from outbreaks, with one having both VEEV complex and MADV RNAs. VEEV complex RNA was found in five suspected dengue cases from disease surveillance. The rRT-PCR assays identified VEEV complex RNA in three Culex (Melanoconion) vomerifer pools, leading to VEEV isolates in two. Phylogenetic analysis revealed the VEEV ID subtype in positive samples. Notably, 11.9% of dengue-like disease patients showed VEEV infections. Together, our rRT-PCR validation in human and mosquito samples suggests that this method can be incorporated into mosquito and human encephalitic alphavirus surveillance programs in endemic regions.

Convergent trends and spatiotemporal patterns of Aedes-borne arboviruses in Mexico and Central America.

BACKGROUND: Aedes-borne arboviruses cause both seasonal epidemics and emerging outbreaks with a significant impact on global health. These viruses share mosquito vector species, often infecting the same host population within overlapping geographic regions. Thus, comparative analyses of the virus evolutionary and epidemiological dynamics across spatial and temporal scales could reveal convergent trends. METHODOLOGY/PRINCIPAL FINDINGS: Focusing on Mexico as a case study, we generated novel chikungunya and dengue (CHIKV, DENV-1 and DENV-2) virus genomes from an epidemiological surveillance-derived historical sample collection, and analysed them together with longitudinally-collected genome and epidemiological data from the Americas. Aedes-borne arboviruses endemically circulating within the country were found to be introduced multiple times from lineages predominantly sampled from the Caribbean and Central America. For CHIKV, at least thirteen introductions were inferred over a year, with six of these leading to persistent transmission chains. For both DENV-1 and DENV-2, at least seven introductions were inferred over a decade. CONCLUSIONS/SIGNIFICANCE: Our results suggest that CHIKV, DENV-1 and DENV-2 in Mexico share evolutionary and epidemiological trajectories. The southwest region of the country was determined to be the most likely location for viral introductions from abroad, with a subsequent spread into the Pacific coast towards the north of Mexico. Virus diffusion patterns observed across the country are likely driven by multiple factors, including mobility linked to human migration from Central towards North America. Considering Mexico's geographic positioning displaying a high human mobility across borders, our results prompt the need to better understand the role of anthropogenic factors in the transmission dynamics of Aedes-borne arboviruses, particularly linked to land-based human migration.

Madariaga and Venezuelan equine encephalitis virus seroprevalence in rodent enzootic hosts in Eastern and Western Panama.

While rodents are primary reservoirs of Venezuelan equine encephalitis virus (VEEV), their role in Madariaga virus (MADV) transmission remains uncertain, particularly given their overlapping geographic distribution. This study explores the interplay of alphavirus prevalence, rodent diversity, and land use within Darien and Western Panama provinces. A total of three locations were selected for rodent sampling in Darien province: Los Pavitos, El Real de Santa Maria and Santa Librada. Two sites were selected in Western Panama province: El Cacao and Cirí Grande. We used plaque reduction neutralization tests to assess MADV and VEEV seroprevalences in 599 rodents of 16 species across five study sites. MADV seroprevalence was observed at higher rates in Los Pavitos (Darien province), 9.0%, 95% CI: 3.6-17.6, while VEEV seroprevalence was elevated in El Cacao (Western Panama province), 27.3%, 95% CI: 16.1-40.9, and El Real de Santa María (Darien province), 20.4%, 95% CI: 12.6-29.7. Species like Oryzomys coesi, 23.1%, 95% CI: 5.0-53.8, and Transandinomys bolivaris, 20.0%, 95% CI: 0.5-71.6 displayed higher MADV seroprevalences than other species, whereas Transandinomys bolivaris, 80.0%, 95% CI: 28.3-99.4, and Proechimys semispinosus, 27.3%, 95% CI: 17.0-39.6, exhibited higher VEEV seroprevalences. Our findings provide support to the notion that rodents are vertebrate reservoirs of MADV and reveal spatial variations in alphavirus seropositivity among rodent species, with different provinces exhibiting distinct rates for MADV and VEEV. Moreover, specific rodent species are linked to unique seroprevalence patterns for these viruses, suggesting that rodent diversity and environmental conditions might play a significant role in shaping alphavirus distribution.

Dispersal patterns and influence of air travel during the global expansion of SARS-CoV-2 variants of concern.

The Alpha, Beta, and Gamma SARS-CoV-2 variants of concern (VOCs) co-circulated globally during 2020 and 2021, fueling waves of infections. They were displaced by Delta during a third wave worldwide in 2021, which, in turn, was displaced by Omicron in late 2021. In this study, we use phylogenetic and phylogeographic methods to reconstruct the dispersal patterns of VOCs worldwide. We find that source-sink dynamics varied substantially by VOC and identify countries that acted as global and regional hubs of dissemination. We demonstrate the declining role of presumed origin countries of VOCs in their global dispersal, estimating that India contributed <15% of Delta exports and South Africa <1%-2% of Omicron dispersal. We estimate that >80 countries had received introductions of Omicron within 100 days of its emergence, associated with accelerated passenger air travel and higher transmissibility. Our study highlights the rapid dispersal of highly transmissible variants, with implications for genomic surveillance along the hierarchical airline network.

Spatiotemporal dynamics and recurrence of chikungunya virus in Brazil: an epidemiological study.

BACKGROUND: Chikungunya virus (CHIKV) is an Aedes mosquito-borne virus that has caused large epidemics linked to acute, chronic, and severe clinical outcomes. Currently, Brazil has the highest number of chikungunya cases in the Americas. We aimed to investigate the spatiotemporal dynamics and recurrence pattern of chikungunya in Brazil since its introduction in 2013. METHODS: In this epidemiological study, we used CHIKV genomic sequencing data, CHIKV vector information, and aggregate clinical data on chikungunya cases from Brazil. The genomic data comprised 241 Brazilian CHIKV genome sequences from GenBank (n=180) and the 2022 CHIKV outbreak in Ceará state (n=61). The vector data (Breteau index and House index) were obtained from the Brazilian Ministry of Health for all 184 municipalities in Ceará state and 116 municipalities in Tocantins state in 2022. Epidemiological data on laboratory-confirmed cases of chikungunya between 2013 and 2022 were obtained from the Brazilian Ministry of Health and Laboratory of Public Health of Ceará. We assessed the spatiotemporal dynamics of chikungunya in Brazil via time series, mapping, age-sex distribution, cumulative case-fatality, linear correlation, logistic regression, and phylogenetic analyses. FINDINGS: Between March 3, 2013, and June 4, 2022, 253 545 laboratory-confirmed chikungunya cases were reported in 3316 (59·5%) of 5570 municipalities, mainly distributed in seven epidemic waves from 2016 to 2022. To date, Ceará in the northeast has been the most affected state, with 77 418 cases during the two largest epidemic waves in 2016 and 2017 and the third wave in 2022. From 2016 to 2022 in Ceará, the odds of being CHIKV-positive were higher in females than in males (odds ratio 0·87, 95% CI 0·85-0·89, p<0·0001), and the cumulative case-fatality ratio was 1·3 deaths per 1000 cases. Chikungunya recurrences in the states of Ceará, Tocantins (recurrence in 2022), and Pernambuco (recurrence in 2021) were limited to municipalities with few or no previously reported cases in the previous epidemic waves. The recurrence of chikungunya in Ceará in 2022 was associated with a new East-Central-South-African lineage. Population density metrics of the main CHIKV vector in Brazil, Aedes aegypti, were not correlated spatially with locations of chikungunya recurrence in Ceará and Tocantins. INTERPRETATION: Spatial heterogeneity of CHIKV spread and population immunity might explain the recurrence pattern of chikungunya in Brazil. These results can be used to inform public health interventions to prevent future chikungunya epidemic waves in urban settings. FUNDING: Global Virus Network, Burroughs Wellcome Fund, Wellcome Trust, US National Institutes of Health, São Paulo Research Foundation, Brazil Ministry of Education, UK Medical Research Council, Brazilian National Council for Scientific and Technological Development, and UK Royal Society. TRANSLATION: For the Portuguese translation of the abstract see Supplementary Materials section.

Dynamics of Early Establishment of SARS-CoV-2 VOC Omicron Lineages in Minas Gerais, Brazil.

Brazil is one of the nations most affected by Coronavirus disease 2019 (COVID-19). The introduction and establishment of new virus variants can be related to an increase in cases and fatalities. The emergence of Omicron, the most modified SARS-CoV-2 variant, caused alarm for the public health of Brazil. In this study, we examined the effects of the Omicron introduction in Minas Gerais (MG), the second-most populous state of Brazil. A total of 430 Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) samples from November 2021 to June 2022 from Belo Horizonte (BH) city were sequenced. These newly sequenced genomes comprise 72% of all previously available SARS-CoV-2 genomes for the city. Evolutionary analysis of novel viral genomes reveals that a great diversity of Omicron sublineages have circulated in BH, a pattern in-keeping with observations across Brazil more generally. Bayesian phylogeographic reconstructions indicate that this diversity is a product of a large number of international and national importations. As observed previously, São Paulo state is shown as a significant hub for viral spread throughout the country, contributing to around 70% of all viral Omicron introductions detected in MG.

Global Expansion of SARS-CoV-2 Variants of Concern: Dispersal Patterns and Influence of Air Travel.

In many regions of the world, the Alpha, Beta and Gamma SARS-CoV-2 Variants of Concern (VOCs) co-circulated during 2020-21 and fueled waves of infections. During 2021, these variants were almost completely displaced by the Delta variant, causing a third wave of infections worldwide. This phenomenon of global viral lineage displacement was observed again in late 2021, when the Omicron variant disseminated globally. In this study, we use phylogenetic and phylogeographic methods to reconstruct the dispersal patterns of SARS-CoV-2 VOCs worldwide. We find that the source-sink dynamics of SARS-CoV-2 varied substantially by VOC, and identify countries that acted as global hubs of variant dissemination, while other countries became regional contributors to the export of specific variants. We demonstrate a declining role of presumed origin countries of VOCs to their global dispersal: we estimate that India contributed <15% of all global exports of Delta to other countries and South Africa <1-2% of all global Omicron exports globally. We further estimate that >80 countries had received introductions of Omicron BA.1 100 days after its inferred date of emergence, compared to just over 25 countries for the Alpha variant. This increased speed of global dissemination was associated with a rebound in air travel volume prior to Omicron emergence in addition to the higher transmissibility of Omicron relative to Alpha. Our study highlights the importance of global and regional hubs in VOC dispersal, and the speed at which highly transmissible variants disseminate through these hubs, even before their detection and characterization through genomic surveillance. Highlights: Global phylogenetic analysis reveals relationship between air travel and speed of dispersal of SARS-CoV-2 variants of concern (VOCs)Omicron VOC spread to 5x more countries within 100 days of its emergence compared to all other VOCsOnward transmission and dissemination of VOCs Delta and Omicron was primarily from secondary hubs rather than initial country of detection during a time of increased global air travelAnalysis highlights highly connected countries identified as major global and regional exporters of VOCs.

Global disparities in SARS-CoV-2 genomic surveillance.

Genomic sequencing is essential to track the evolution and spread of SARS-CoV-2, optimize molecular tests, treatments, vaccines, and guide public health responses. To investigate the global SARS-CoV-2 genomic surveillance, we used sequences shared via GISAID to estimate the impact of sequencing intensity and turnaround times on variant detection in 189 countries. In the first two years of the pandemic, 78% of high-income countries sequenced >0.5% of their COVID-19 cases, while 42% of low- and middle-income countries reached that mark. Around 25% of the genomes from high income countries were submitted within 21 days, a pattern observed in 5% of the genomes from low- and middle-income countries. We found that sequencing around 0.5% of the cases, with a turnaround time <21 days, could provide a benchmark for SARS-CoV-2 genomic surveillance. Socioeconomic inequalities undermine the global pandemic preparedness, and efforts must be made to support low- and middle-income countries improve their local sequencing capacity.

First Report of Wenzhou sobemo-like virus 4 in Aedes albopictus (Diptera: Culicidae) in Latin America.

Insect-specific viruses (ISVs) are viruses that replicate exclusively in arthropod cells. Many ISVs have been studied in mosquitoes as many of them act as vectors for human etiological agents, such as arboviruses. Aedes (Stegomyia) albopictus is an important potential vector of several arboviruses in Brazil, such as dengue (DENV), Zika (ZIKV) and chikungunya (CHIKV). The development of next-generation sequencing metagenomics has enabled the discovery and characterization of new ISVs. Ae. albopictus eggs were collected using oviposition traps placed in two urban parks in the city of São Paulo, Brazil. The Aedes albopictus females were divided into pools and the genetic material was extracted and processed for sequencing by metagenomics. Complete genomes of ISV Wenzhou sobemo-like virus 4 (WSLV4) were obtained in three of the four pools tested. This is the first detection of ISV WSLV4 in Ae. albopictus females in Latin America. Further studies on ISVs in Ae. albopictus are needed to better understand the role of this species in the dynamics of arbovirus transmission in the Americas.

Predicting SARS-CoV-2 Variant Spread in a Completely Seropositive Population Using Semi-Quantitative Antibody Measurements in Blood Donors.

SARS-CoV-2 serologic surveys estimate the proportion of the population with antibodies against historical variants, which nears 100% in many settings. New approaches are required to fully exploit serosurvey data. Using a SARS-CoV-2 anti-Spike (S) protein chemiluminescent microparticle assay, we attained a semi-quantitative measurement of population IgG titers in serial cross-sectional monthly samples of blood donations across seven Brazilian state capitals (March 2021−November 2021). Using an ecological analysis, we assessed the contributions of prior attack rate and vaccination to antibody titer. We compared anti-S titer across the seven cities during the growth phase of the Delta variant and used this to predict the resulting age-standardized incidence of severe COVID-19 cases. We tested ~780 samples per month, per location. Seroprevalence rose to >95% across all seven capitals by November 2021. Driven by vaccination, mean antibody titer increased 16-fold over the study, with the greatest increases occurring in cities with the highest prior attack rates. Mean anti-S IgG was strongly correlated (adjusted R2 = 0.89) with the number of severe cases caused by Delta. Semi-quantitative anti-S antibody titers are informative about prior exposure and vaccination coverage and may also indicate the potential impact of future SARS-CoV-2 variants.

SARS-CoV-2 antibody dynamics in blood donors and COVID-19 epidemiology in eight Brazilian state capitals: A serial cross-sectional study.

Background: The COVID-19 situation in Brazil is complex due to large differences in the shape and size of regional epidemics. Understanding these patterns is crucial to understand future outbreaks of SARS-CoV-2 or other respiratory pathogens in the country. Methods: We tested 97,950 blood donation samples for IgG antibodies from March 2020 to March 2021 in 8 of Brazil's most populous cities. Residential postal codes were used to obtain representative samples. Weekly age- and sex-specific seroprevalence were estimated by correcting the crude seroprevalence by test sensitivity, specificity, and antibody waning. Results: The inferred attack rate of SARS-CoV-2 in December 2020, before the Gamma variant of concern (VOC) was dominant, ranged from 19.3% (95% credible interval [CrI] 17.5-21.2%) in Curitiba to 75.0% (95% CrI 70.8-80.3%) in Manaus. Seroprevalence was consistently smaller in women and donors older than 55 years. The age-specific infection fatality rate (IFR) differed between cities and consistently increased with age. The infection hospitalisation rate increased significantly during the Gamma-dominated second wave in Manaus, suggesting increased morbidity of the Gamma VOC compared to previous variants circulating in Manaus. The higher disease penetrance associated with the health system's collapse increased the overall IFR by a minimum factor of 2.91 (95% CrI 2.43-3.53). Conclusions: These results highlight the utility of blood donor serosurveillance to track epidemic maturity and demonstrate demographic and spatial heterogeneity in SARS-CoV-2 spread. Funding: This work was supported by Itaú Unibanco 'Todos pela Saude' program; FAPESP (grants 18/14389-0, 2019/21585-0); Wellcome Trust and Royal Society Sir Henry Dale Fellowship 204311/Z/16/Z; the Gates Foundation (INV- 034540 and INV-034652); REDS-IV-P (grant HHSN268201100007I); the UK Medical Research Council (MR/S0195/1, MR/V038109/1); CAPES; CNPq (304714/2018-6); Fundação Faculdade de Medicina; Programa Inova Fiocruz-CE/Funcap - Edital 01/2020 Number: FIO-0167-00065.01.00/20 SPU N°06531047/2020; JBS - Fazer o bem faz bem.

Using multiple sampling strategies to estimate SARS-CoV-2 epidemiological parameters from genomic sequencing data.

The choice of viral sequences used in genetic and epidemiological analysis is important as it can induce biases that detract from the value of these rich datasets. This raises questions about how a set of sequences should be chosen for analysis. We provide insights on these largely understudied problems using SARS-CoV-2 genomic sequences from Hong Kong, China, and the Amazonas State, Brazil. We consider multiple sampling schemes which were used to estimate Rt and rt as well as related R0 and date of origin parameters. We find that both Rt and rt are sensitive to changes in sampling whilst R0 and the date of origin are relatively robust. Moreover, we find that analysis using unsampled datasets result in the most biased Rt and rt estimates for both our Hong Kong and Amazonas case studies. We highlight that sampling strategy choices may be an influential yet neglected component of sequencing analysis pipelines.

Tracking the turnover of SARS-CoV-2 VOCs Gamma to Delta in a Brazilian state (Minas Gerais) with a high-vaccination status.

The emergence and global dissemination of Severe Acute Respiratory Syndrome virus 2 (SARS-CoV-2) variants of concern (VOCs) have been described as the main factor driving the Coronavirus Disease 2019 pandemic. In Brazil, the Gamma variant dominated the epidemiological scenario during the first period of 2021. Many Brazilian regions detected the Delta variant after its first description and documented its spread. To monitor the introduction and spread of VOC Delta, we performed Polymerase Chain Reaction (PCR) genotyping and genome sequencing in ten regional sentinel units from June to October 2021 in the State of Minas Gerais (MG). We documented the introduction and spread of Delta, comprising 70 per cent of the cases 8 weeks later. Comparing the viral loads of the Gamma and Delta dominance periods, we provide additional evidence that the latter is more transmissible. The spread and dominance of Delta did not culminate in the increase in cases and deaths, suggesting that the vaccination may have restrained the epidemic growth. Analysis of 224 novel Delta genomes revealed that Rio de Janeiro state was the primary source for disseminating this variant in the state of MG. We present the establishment of Delta, providing evidence of its enhanced transmissibility and showing that this variant shift did not aggravate the epidemiological scenario in a high immunity setting.

Context-specific emergence and growth of the SARS-CoV-2 Delta variant.

The SARS-CoV-2 Delta (Pango lineage B.1.617.2) variant of concern spread globally, causing resurgences of COVID-19 worldwide1,2. The emergence of the Delta variant in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions. Here we analyse 52,992 SARS-CoV-2 genomes from England together with 93,649 genomes from the rest of the world to reconstruct the emergence of Delta and quantify its introduction to and regional dissemination across England in the context of changing travel and social restrictions. Using analysis of human movement, contact tracing and virus genomic data, we find that the geographic focus of the expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced more than 1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers reduced onward transmission from importations; however, the transmission chains that later dominated the Delta wave in England were seeded before travel restrictions were introduced. Increasing inter-regional travel within England drove the nationwide dissemination of Delta, with some cities receiving more than 2,000 observable lineage introductions from elsewhere. Subsequently, increased levels of local population mixing-and not the number of importations-were associated with the faster relative spread of Delta. The invasion dynamics of Delta depended on spatial heterogeneity in contact patterns, and our findings will inform optimal spatial interventions to reduce the transmission of current and future variants of concern, such as Omicron (Pango lineage B.1.1.529).